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Abstract 

We present an evaluation of the strong couplings J D'^*'> D'^*'^ and J D^*'> D'-*'> t: by an effective 
field theory of quarks and mesons. These couplings are necessary to calculate n + J/ip 
£)(*) _(_ jji*) cross sections, an important background to the J/^ suppression signal in the 
quark-gluon plasma. We write down the general effective Lagrangian and compute the 
relevant couplings in the soft pion limit and beyond. 
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We present an evaluation of the strong couplings JD ^' D *' and JD *' D *'-k by an effective field 
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This paper is devoted to the study of the strong couphngs of J low mass charmed mesons and pions. The 
interest of this study stems from the possibility that J /ijj absorption processes of the following type 



play an important role in the relativistic heavy ion scattering. Since a decrease of the J production might signal the 
formation of Quark Gluon Plasma (QGP) in a heavy ion collision, it is useful to have reliable estimates of the cross 
sections for the processes that provide an alternative way to reduce the J/ip production rate. Previous studies of 
these effects can be found in il|,,^,3||. The relevant couphngs needed to compute ^ are depicted in Fig.Q] 
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I. INTRODUCTION 
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FIG. 1: Couplings involved in the tree level calculations of the processes ITl . 



Besides the DD*it coupling;, see Fig.^, whose couphng constant go* Dtt has been theoreticahy estimated Pj, and 
experimentally investigated to compute the amplitudes ^ one would need also the JD'^*^D''*\ see Fig.^, and the 
J D^*") D'^*'> TT couplings, see Fig.^. In an effective Lagrangian approach the latter couplings provide direct four-body 
interactions, while the former enter the amplitude via tree diagrams with the exchange of a charmed particle in 
the t— channel. These couplings have been estimated by different methods, that are, in our opinion, unsatisfactory. 
For example the use of the SU4 symmetry puts on the same footing the heavy quark c and the light quarks, which is at 
odds with the results obtained within the Heavy Quark Effective Theory (HQET) , where the opposite approximation 
rric ^ ^QCD is used (for a short review of HQET see 0)- Similarly, the rather common approach based on the 
Vector Meson Dominance (VMD) should be considered critically, given the large extrapolation — "^J/V' ^^^^ 
is involved. A different evaluation, based on QCD Sum Rules can be found in and presents the typical theoretical 
uncertainties of this method. In this note we will use a different approach, based on the Constituent Quark Meson 
(CQM) model 8], 9], which takes into account explicitly the HQET symmetries. We gave a preliminary report of the 
present study in where we presented the results for some of the couplings Here we complete the analysis and 
compare our results with the existing literature. 



II. CONSTITUENT QUARK MESON MODEL 

The CQM model is a quark-meson model arising from an extension of Ref. ^3 to exploit fully HQET and chiral 
symmetries in the interactions of heavy and light mesons. A survey of these methods can be found in 0,0 ■ Here 
we limit to those aspects of the model that are relevant for the interactions of J/ip, low mass charmed mesons and 
pions. The model is an effective field theory whose effective fields are light and heavy quarks as well as light and 
heavy mesons. The Feynman rules for the model are explicitly written down in 0, 0. The transition amplitudes 
containing light /heavy mesons in the initial and final states as well as the couplings of the heavy mesons to hadronic 
currents are computable via quark loop diagrams where mesons enter as external legs. The model is relativistic and 
incorporates, besides the heavy quark symmetries, also the chiral symmetry of the light quark sector. 

To show an example of the calculation in the CQM model we consider the Isgur-Wise (IW) function defined e.g. by 

{Dip')\crc\D{p)) = mDa^){v + v'Y (2) 

where = muv^, p' ^ = rriDv' ^ and lu = v ■ v' . We note that the Isgur-Wise function obeys the normalization 
condition ^(1) = 1, arising from the flavor symmetry of the HQET. The explicit definition of the Isgur-Wise form 
factor is: 

{H{v')\cj^c\H{v)) = -e(c^)Tr {H^,,H) . (3) 
Here H is the multiplet containing both the D and the D* mesons 0: 

^ = ^^^^(-^75+7-^') , (4) 
and P5, are annihilation operators for the charmed mesons. One gets 

u , 7 3t f 4 Tr [(7-^-1- m)75(l + 7 ■ v'^^^jl + 7 ■ ^7)75] , , 

where v and v' are the 4— velocities of the two heavy quarks that are equal, in the infinite quark mass limit, to the 
hadron velocities, 



1 -|- 7 • 11 
2(z;-£-f Ah) 



(6) 



is the heavy quark propagator of the HQET, and A// — mo — riic — v ■ k m the limit rric -^00; fc is the meson residual 
momentum, defined by p^ — nicV^^ + . The numerical value of A// is in the range 0.3 — 0.5 GeV 0. If we consider 
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a D* meson instead of a Z?, a factor —75 must be substituted by 7 • e, e being the polarization of D*. The constant 
Zh is the heavy meson field wavefunction renormalization constant giving the strength of the quark-meson coupling 
(more precisely the coupling is ^/Zhitiu). Zh is computed and tabulated in 8]. One gets for the IW function: 

2 



Ciio) = Zh 
where the li integrals are given by: 

^3 (A) = 

/5(Ai,A2,w) = 



l+U! 



iNr 



(fc2 - rn?){v • fc + A + ze) 



167r4 

_N, 
16 

167r4 
1 

dx 



rcg 



2 A 

m+— ^/5(Aj^,Ah,w) 
1 + w 



(7) 



(8) 



- m?){v • fc + Ai + ie)(i;' • fc + A2 + ie) 
1 







l + 2x^{l-uj) + 2x{uj-l) 
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16^3/2 7^/^ 
6 



ds cr e-'*('"'-'^') s-1/2 (1 + erf (dVi)) + 



167r2 



1/A2 



where 



CT(a:, Ai, A2,a;) 



Ai (l-a;) + A2a; 



^1 + 2 [uo -I) X + 2 [l-uj) 



(9) 



(10) 



The ultraviolet cutoff A, the infrared cutoff /i and the light constituent mass m are fixed in the model [3 to be 
A = 1.25 GeV, /i = 0.3 GeV and m = 0.3 GeV. Other integrals to be used later are 



h 

/4(A) 



167r4 
167r4 

Nr. 



Nr. 



r 0, 



(fc^ — m^)^ 16 TT^ 

= d^ 

(fc2 - ni?Y{v ■ k + A + ie) 

1/^^ d,s 



2 2 
m TO 

A^V 



-s(m^-A^) 



The calculation of the qd'Ot! coupling constant for the matrix element 

{iT+{q)D^{p)\D*+{p',e))=zgD.D.e-q 



[1 + erf (A\/i)] 



(11) 



(12) 



(13) 



in the CQM proceeds along similar lines and can be found in ^8]. We reproduce for reference's sake since it is an 
important element of the amplitudes In the soft pion limit (spl) one gets |8|: 



5^^i^^ = 13±l 



(14) 



The experimental situation is as follows: CLEO results Q give g — 0.59±0.07±0.01 where g is related to the constant 
in eq. ((TH) by g'^j^l^,^ = 2mDg/ fir (/tt ~ 130 MeV); on the other hand, in general, QCD sum rules predict smaller 
values, see for a review 



4 



III. JD^*^D^*^ COUPLINGS 

The calculation of the Isgur-Wise we have described above is a crucial ingredient to the computation of the J D^*'> D'^*'> 
vertexes of Fig. It corresponds to the evaluation of the l.h.s. of Fig. |21 while, via a VMD ansatz, the r.h.s. gives 
the desired coupling. Concerning the use of the VMD for the charm system one has to note that it is not based on the 
hypothesis that all higher order resonances give contributions smaller than the 3/ip^ but on the fact that the higher 
states give contributions of alternating sign that tend to cancel. This sign difference follows from an evaluation via 
the saddle point method of the WKB approximation • 




J ^ D 



D 



FIG. 2: The Vector Meson Dominance equation giving the coupling of J/ip with D,D* in terms of the Isgur-Wise function The 
function J on the l.h.s. is computed by a diagram with a quark loop. The coupling of each D^*^ meson to quarks is given by y/Z~[pmo. 

The Isgur-Wise function can be computed for any value of and not only in the region a; > 1, which is experimentally 
accessible via the semileptonic B — s- D^*^ decays; lo is related to the meson momenta by 



(15) 



where Pi, P2 = momenta of the two D resonances. 

Let us now consider the r.h.s. of the equation depicted in Fig.|21 For the coupling of J/ip to the current we use the 
matrix element 



(0|c7^c| J(g, 77)) =/jmje^ 



(16) 



with /./ = 0.405 ± 0.014 GeV. As to the strong couplings J£)(*)L>(*), the model in Fig. H gives the following effective 
Lagrangians 



CjDD = ig.jDD {odi^D^ .r , 



(17) 
(18) 



JD*D* 



(19) 



Here !)(*)£)(*) can be any of the pairs D°^*^D°^*\ £)+(*)£)-(*) or oi*'' oi*^ (neglecting SU3 breaking effects). As a 
consequence of the spin symmetry of the HQET we find: 



gjD-D' — gjDD ■ 



gjDD' 



gjDD 

mo 



while the VMD ansatz gives: 



2 2 

gjDDiPi, P2^P )= jj^j ^^^1 ■ 



(20) 



(21) 
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Since gjoD has no zeros, eq. (|21ll shows that ^ has a pole at = m?j, which is what one expects on the basis of 
dispersion relations arguments. The CQM evaluation of ^ does show a strong peak for w (2toc)^, even though, due 
to O (1/toc) effects, the location of the singularity is not exactly at = w?j. 




FIG. 3: The dependence of <; = gjDoi^^t P^)' showing the almost complete cancellation between the pole of the Isgur-Wise 
function and the kinematical zero. Units are GeV^ for p^. 

This is shown in Fig. |3 where we plot gjooiPi, pi: -P^) ^or on shell D mesons, as a function of (the plot is 
obtained for Ah = 0.4 GeV and Zh = 2.36 GeV~^). For in the range (0,4) GeV^, gjoD is almost flat, with a 
value 

=8.0±0.5 . (22) 

For larger values of p^ the method is unreliable due to the above-mentioned incomplete cancellation between the 
kinematical zero and the pole. Therefore, we extrapolate the smooth behavior of gjoD in the small p^ region up to 
p^ = m?j and assume the validity of the result H22() also for on-shell J/ij] mesons. On the other hand in the Pi,;^ 
variables we find a behavior compatible with that produced by a smooth form factor. Also the result from Ref. [ij 
in the Table, is based on a VMD assumption; a previous determination based on the same assumption is p!^ . with 
similar results (gjoD — 7.7). 

In Table 1 we compare our results with those of other authors. We observe that our results for gjDD a-nd gjDD* 
agree with the outcomes of Ref. 13] and with the QCD sum rule analysis of 0; in particular the smooth behavior of 
the form factor found in for gjoD agrees with our result. This is not surprising, as [31 uses a VMD model as well. 

As for the QCD sum rules calculation, it involves a perturbative part and a non perturbative contribution, which 
is however suppressed. The perturbative term has its counterpart in CQM in the loop calculation of Fig. Eland the 
overall normalization should agree as a consequence of the Luke's theorem. On the other hand we differ from Ref. 
for a sign and from both in sign and in magnitude. The sign difference may be due to an overall phase in the 
J/^E" wavefunction. It is however of no effect in the computation of the J/'^ absorption cross section, which is the 
main application of the present calculation. Finally Ref. 13] obtains a value for gjDD' using results from the decay 
J/f/; pn. This seems to us too strong an assumption due to the fact that J/^p — > pn could proceed via gluonic 
decay of the J/ip, which is not the case for J/i/; — > DD* . 



Coupling 


Our work 


Ref. \l] 


Ref. \IA] 


Ref. [13] 


Ref. [2 


gjDD 


8.0 ±0.5 


-7.64 


-4.93 


7.71 


7.36 


gjDD' (GeY-') 


4.05 ±0.25 






8.02 ±0.62 




gjD'D* 


8.0 ±0.5 


-7.64 


-4.93 


7.71 





Table 1. Comparison of theoretical results for the couplings gjDD, gjDD* and gjo'D'- Ref. Q and Ref. use a VMD model similar 
to the one used in the present paper for the couplings gjoD^ gjD*D* ■ For gjoD* Ref. Il3l uses VMD together with data from relativistic 
quark model |15| to get the coupling of a hadronic current to D and D* . Ref. |14| uses a chiral model to compute the coupling constants 
gjDD and gjD*D* ■ The coupling gjDD' is not included in Ref. Ref. is based on QCD sum rules; the result we report for gjDD is 
computed at the same value = 2 GeV^ as in our work. 
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IV. JD^*^D^*K COUPLINGS 

Let us now consider the J D^*") D^*^ t: couplings of Fig. We write the effective Lagrangians for the Jn^D^D^ 
coupHng (other couphngs can be obtained by use of symmetries): 

CjDD. = -^e^''"''j^{d,7T)do.Ddi3D, (23) 
mo 

^JD'D. = ~Md,.n){g,[g'^''Dld^D + 2Dd''D*'']~g2[D*''d'' + Dd''D*''] 

+ -^[{d''D*^)d^'dxD-id''dxD)d^'D*^]\ + h.c. , (24) 
irijj J 

^ 771 J-) 



97 



m 



D 



m 



D 



- ^ [epo/3A(a"a.i?*'')9^a,5*'' + e„p^A(a"9.i?*'')a^9p5*'3] 

+ ^ [epa/3A (a"i?*o ap9,a^^*^' + e^p^^ (a"/?*") a^a^a,^*'^] | . (25) 

While in these formulae 13 coupling constants appear, the number of the independent couplings is only 5. As a matter 
of fact they can be written in terms of the independent couplings A\^Ai^A/^^B^K defined by the formulae 



167r4 



3i 
167r4 
3i 



(^2- 




hg)2 


— m?] [v 
+ A2 


^ + A)(^;' • 


^ + A) 

) +A4 




m2)[(£H 


h9)2 


— to2] (^y 
1 


£ + A)(t;' • 


^ + A) 


(^2- 


m2)[(£H 




— m2](u 


i + A){v' ■ 


^ + A) 



167r4 

The dependence of the gk on these couplings is in Table 2 



d'^f - ^ B iv" + v'f^) 
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Coupling 


CQM formula 


On 




P[{4:m + 2AH)B -m^K] 




y i 


(3 [2Ai + 2A2 + 2ujA4 - AmB + rr 




g2 


13 [2Ai - 


f 2ujA2 + 2A4, - 2(1 + Lo)mB 




ifo 




2(3 [A2 - ^4 - mB] 




94 


f3[2A, 


+ 2A2 + 2uoAi - (1 + uj)mB - 




95 




/3 [2A2 - mB] 




96 




(3 [2Ai + 2Ai ^mB + m^K] 




97 




l3mB 




9'7 




51 + 257 




5? 




51 + 97 




58 




2(3[{l-uj)A2- Ai] 




59 




2[3Ai 




9'9 




95 + 97 





Table 2. Relations between the coupling constants g^. and the constants Ai, A2, A4, B and K. Units are GeV ^; ui is defined in ea. ll5i . 



vith 



-Pi- 



and = 2 GeV^ for go i 9i 1 94 1 97 and p-^ = 5 GcV^ for 92 , 93 , 95 , 



,99- 



In this table 



(3 



P 



fyrfjmj 



Zh ; 



(27) 



and explicit formulae for Ai , A2 , A4 , B , K can be found in Section These results have been obtained by a VMD 
ansatz similar to Fig. |21 but now the l.h.s is modified by the insertion of a soft pion on the light quark line (with a 
coupling 9^ //tt 7/^75)- Let us discuss in some detail one of these couplings, go. The numerical results for on-shell D 
mesons, in the soft pion limit as a function of the J/^p virtuality show a behavior similar to that of Fig. |21 By the 
same arguments used to determine Qjdd in Fig.|31we choose — 2 GeV^ and we get 



5o 
ruD 



125 ± 15 GeV" 



(soft pion limit) . 

In order to include hard pion effects we write the general formula 

9k{\q.\)=9k{0)fk{\q.\), 



(28) 



(29) 



where fk{\(lTv\) is a form factor. We will discuss it in the next section. For the time being we report the values of all 
the coupling constants in the soft pion limit in Table 3. 



Coupling 


Results (GcV-2) 


Coupling 


Resuhs (GeV-2) 


5o(0) 


+234 ± 45 


55(0) 


-274 ± 49 


51 (0) 


-235 ± 38 


56(0) 


+ 104 ± 16 


52(0) 


-126 ± 19 


57(0) 


+30 + 5 


53(0) 


-252 ± 38 


58(0) 


-106 + 22 


54(0) 


-165 ±25 


59(0) 


-63 + 37 



Table 3. Results for the coupling constants 9j(|9ir|) in the soft pion limit 



In computing this table we have adopted a criterion of stability in analogous to the one used for gQ. For 
5o J 51 , 54 , 57 we find stability at the J/ ip virtuality p^ = 2 GeV^ . For 52 , 53 , 55 , 58 > 59 we find stability at = 5 
GeV^. The technical reason for this difference is that, in the latter case, the equations do not determine the five 
constants for p^ w 0; therefore the stability region lies around the center of the p^ interval (0,TOj). For g^ we do not 
find stability and we derive it by consistency equations derived from Table 2. To the error arising from the stability 
analysis we have added a further theoretical uncertainty of +15% in quadrature. 
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An attempt to compute quadrilinear couplings via SU4 symmetry relations can be found in [l| and in p^ . For 
example, the result for go obtained in is 30 GeV~^. The difference with Table 3 is due to the large SU4 violations 
of our model {nic >> Tnu,ms). 



V. FORM FACTORS 



The coupling constants gj can be expressed in terms of the constants Aj, B and K using the results of Table 2. 
These constants, for ^ 0, are expressed in terms of parametric integrals Ij as follows: 



A2 
Ai 



Cljl - LU^) - 2C3 + 2C2UJ 

2(1 -t^2) ' 

Cijiu^-l)- 6C2CJ + 2C3(2 + Lj^) 
2{uj^ - 1)2 ' 



where 



2(cj2 - 1)2 

1 

B{q.) = -— dx[h{A'{x,q^))- A{x,q^)Kix,q^)] , (30) 
1 Jo 

Kiq^) = / dxK{x,q^), (31) 



1 

Ciiq^) = I dx[l5{A{x,q^),A'{x,q^),u;)+m'^K{x,q^)] 



C2(g7r) = -I2- / dx[A'{x,q.^)l4{A'{x,q„)) + A{x,q.^)l4{A{x,q.^)) 

JQ 

- A{x,q^)A'{x,q^)K{x,q^)] , 
C^^q-n) = -L0I2- I dx[{A{x,q^) + u!A'{x,qT,))h{A'{x,q^))A^{x,q^)K{x,q^)] 



K{x,q^) - ^/5(A(a;,(7,),A'(x,g^),w), (32) 
and we have defined 

A(a;,9^) = A - g^x , (33) 

A'{x,q^) = A-q^xu . (34) 

To compute the integrals we have applied the Feynman trick with the shift £ + gx — > £' where the pion momentum 
is q^ = {q-n,Q, 0, g^). uj is computed by eq. p5|) with the appropriate value of p^^ according to the discussion above. 
Moreover we have used the approximation: 

w'^^ojw", (35) 

which has the correct normalization at w = 1 and also satisfies the constraint oj — v -v' . Using H35() we assume that at 
least one of the two charmed resonances is off-shell, simplifying considerably the numerical computation. A numerical 
calculation of the integrals in these equations leads to a general fit of the form factors as follows: 

hm) = — \^ , (36) 
1 + ^ 

TOfe 

with approximately the same value for all the form factors: 

nik =0.20 ±0.05 GeV . (37) 
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g 




FIG. 4: The dependence of g(q) = go* D-n(<l) / dDDiri'^) on the pion momentum q. 



It is useful to compute the corrections to the soft pion limit also for the DD*tt coupling constants whose value 
for — !■ was computed in ^sj and reported in eq. (|14|l . Using the same technique employed above we make the 
substitution 



1 



1 



(38) 



and compute the integrals appearing in eq. (63) of [g] by the Feynman method as above. The result of this calculation 
is plotted in Fig. 0] The dependence can be fitted by a formula similar to eq. (|36|l with a mass = 0.37 GeV: 



go'D-jiilq-Kl) = 



spl 
9d'Dtt 

1 + P 



(39) 



It is useful to compare this result with that found in Ref. for the same quantity computed in a QCD inspired 
potential model. In that case one finds again a similar form factor with = 0.20 ±0.10 GeV. These two determina- 
tions are sufficiently compatible and induce us to be confident that the method we have used to get the results H3t)|) 
and II37II is reliable. It can be noted that considering a different extrapolation, i.e. in the pion virtuality Q^, as in 
Ref. [131, a smoother behavior is obtained, which should be relevant in the calculation of the J/^P absorption cross 
section. 



VI. CONCLUSIONS 

As discussed in 0, but see also 0, the leading contributions to the current matrix element {H{v')TT\c'y'^c\H{v)) in 
the soft pion limit (spl) are the pole diagrams. The technical reason is that, in the spl, the reducing action of a pion 
derivative in the matrix element is compensated in the polar diagrams by the effect of the denominator that vanishes 
in the combined limit q,^ 0, rric — + 00. Let us now compare this result with the effective JDDtt coupling obtained 
by a polar diagram with an intermediate D* state. We get in this case 

polar _ gjDD'gD'DTrmD , . 

% ~ ^ ■ (40) 

2(77r • PD 

This expression dominates over the result H28|l for pion momenta smaller than 100 MeV. If one restricts the model to 
the soft pion limit (|(f^| < 100 MeV), in spite of the rather large value of go{0) the diagrams containing this coupling 
are suppressed, and one can expect a similar result also for the other couplings. However, to allow the production 
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pair by processes one has to go beyond the spl, since the threshold for the charmed meson pair is 
= 700— 1000 MeV. Our results show that in the CQM model this is indeed possible, by including computed form 
factors as in and l|39|) . Similar form factors were considered in Jj, with a different motivation. Here we have 
shown that the CQM model not only allows their computation, but also gives the general expression for the trilinear 
and quartic couplings of J/ip to charmed mesons and pions. In spite of its model dependent character this seems to 
us an interesting result. Applications to the problem of the J/ip absorption in a nuclear medium will be considered 
elsewhere. 
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